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ABSTRACT 

A new momentum-integral method has been devised for describing the incompress- 
ible flow of a conducting fluid in the entrance of an MHD channel. The method is distinc- 
tive in that an "edge" stress is permitted to exist at the point of intersection of the 
boundary layer and the free stream. A prudent choice of the edge stress variation 
forces the governing differential equation to yield an asymptotic solution with the correct 
wall stress. Accurate analytical solutions are obtained for laminar flow with and with- 
out a magnetic field. Nonmagnetic turbulent flow is also well described by this method. 
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A NEW MOMENTUM- INTEGRAL METHOD FOR TREATING MAGNETOHY- 
DRODYNAMIC AND SIMPLE HYDRODYNAMIC ENTRANCE FLOWS 


by Richard E. Schwirian 
Lewis Research Center 

SUMMARY 

A new momentum-integral method has been devised for describing the incompress- 
ible flow of a conducting fluid in the entrance of a magnetohydrodynamic (MHD) channel. 
The method is distinctive in that an "edge” stress is permitted to exit at the point of in- 
tersection of the boundary layer and the free stream. The edge stress exerts a drag 
force on the free stream and mathematically accounts for the fact that the "free” stream 
becomes less and less free as the fluid enters the channel. A prudent choice of the edge 
stress variation forces the governing differential equation to yield an asymptotic solution 
with the correct wall stress. The edge stress method allows accurate analytical solu- 
tions of the laminar flow momentum integral equations to be obtained for most reason- 
able initial conditions. Flows with and without magnetic fields can be treated. Nonmag- 
netic turbulent flow is also well described by this method. 


INTRODUCTION 

For the designer of liquid-metal magnetohdryodynamic (MHD) generators or elec- 
tromagnetic pumps an accurate knowledge of wall shear stresses is imperative. Indeed, 
the practicality of the MHD generator concepts considered by Elliott (ref. 1) and Patrick 
and Lee (ref. 2) depends largely on keeping the viscous pressure drops across the de- 
vice to a minimum. The most obvious way of meeting this need is to keep such devices 
as short as possible to minimize the time the fluid spends in the high velocity - high 
shear stress power generation section. It is quite possible, therefore, that fully de- 
veloped flow may never be attained in a practical generator. If such is the case, the 
entrance flow problem is of prime importance. 

A number of authors have considered the channel entrance flow problem. Dix 
(ref. 3), Rossow (ref. 4), and Moffatt (ref. 5) consider the laminar flow of a fluid over 



a flat plate perpendicular to which a magnetic field is impressed. Dix’s (ref. 3) solu- 
tions are numerical, Rossow (ref. 4) takes a Blasius series approach, and Moffatt 
(ref. 5) uses a momentum integral method. Moffatt (ref. 5) also considers turbulent 
flow. These solutions are adequate descriptions of entrance flows only when the channel 
half -width a (see fig. 1) is very large relative to the boundary layer thickness 6 and the 
free stream is not significantly accelerated. Bodoia and Osterle (ref. 6), Schlicting 
(ref. 7), Shohet (ref. 8), Roidt and Cess (ref. 9), Moffatt (ref. 10), Maciulaitis and 
Loeffler (ref. 11), and Dhanak (ref. 12) consider the more realistic laminar entrance 



flow problem in which mass conservation demands that the free stream velocity change 
when the boundary layer thickness changes. Bodoia and Osterle (ref. 6) and Schlicting 
(ref. 7) consider the nonmagnetic case with a uniform velocity profile at the entrance. 
Bodoia and Osterle (ref. 6) present a finite -difference technique; Schlicting (ref. 7) 
matches an upstream Blasius series solution to a downstream series solution whose 
zeroth order term is the fully developed profile. Shohet (ref. 8) generalizes the work of 
Bodoia and Osterle (ref. 6) to include magnetic field effects. Similarly, Riodt and Cess 
(ref. 9) generalize Schlicting’ s method (ref. 7) to include magnetic field effects. Integral 
techniques are presented by Moffatt (ref. 10), Maciulaitis and Loeffler (ref. 11), and 
Dhanak (ref. 12). Moffatt (ref. 10) and Maciulaitis and Loeffler (ref. 11) assume simi- 
lar profiles: Moffatt uses a ’’Hartmann” profile; Maciulaitis and Loeffler assume a 
parabolic profile. Dhanak (ref. 12) uses a fourth power polynomial to represent the ve- 
locity profile, applying Karman-Polhausen type conditions to relate the arbitrary con- 
stants. Finally, Maciulaitis and Loeffler (ref. 11) consider turbulent flow, using an in- 
tegral method and a one-seventh power velocity profile. They also consider the case of 
nonmagnetically fully developed entrance conditions. 
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With the exception of the numerical techniques (refs. 3, 6, and 8), all of the pre- 
vious methods have restricted spatial regions of validity. In most cases this means that 
the solutions obtained are valid only for points sufficiently close to the entrance or suf- 
ficiently far downstream. The purpose herein is to introduce an integral technique 
which (1) is valid throughout the channel, (2) gives accurate results with and without a 
magnetic field, (3) can be used with any reasonable entrance profile, including uniform 
entry and nonmagnetically fully developed entry, and (4) is good for laminar and re- 
stricted regions of turbulent flows. All this is made possible by the simple device of 
allowing a viscous edge stress term to exist in the free stream momentum equation. 

This term represents the drag that the boundary layer exerts on the free stream. If the 
boundary layer - free stream model of the flow is to be a valid approximation every- 
where, such a stress must exist to satisfy the end conditions. If it is not included (as 
in refs. 4 to 7 and 9 to 12), the validity of the pressure distributions obtained is confined 
to a small region near the entrance. 

The present investigation is concerned with determining the variation of the edge 
stress and obtaining analytical solutions for the friction factor and pressure defect. 

Both laminar and turbulent flows are considered. 


SYMBOLS 


a 

channel half-width (Y-direction) 


B 

magnetic field strength 


C f 

local friction factor 


<c f > 

average friction factor 


C fB 

Blasius friction factor, 0. 664(ju/pUX) 

1/2 

E 

electric field strength 


f 

"edge" stress function 


K 

load parameter, E/ ( U) B 


m h 

Hartmann number, Ba (a/p.)^^ 


& 

pressure defect, (p - Pq)/p(U)^ 

“lAeH 1 - K > 


hydromagnetic pressure defect, & + ( 

Re 

Reynolds number, 4p(U)a/p 


Re m 

magnetic Reynolds number, p Q 0 'a(U) 

|k — 1 1 
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u 

free stream velocity 

<U> 

average velocity 

u 

velocity in boundary layer (X-direction) 

u* 

dimensionless velocity, u/<U) 

V 

transverse velocity (Y-direction) 

X 

coordinate parallel to flow 

X e 

entrance length 

X 

X/a 

X 1 

4 x/Re 

x 2 

4 M|x/Re 

Y 

coordinate parallel to magnetic field 

y 

Y/a 

z 

coordinate perpendicular to X and Y 

at a ’ 

L’ L 

see eq. (21) 

a T 

see eq. (34) 

6 

boundary layer thickness 

6* 

displacement thickness 

e 

dimensionless boundary layer thickness, 

e* 

e / £ oc 

V 

Y/6 

e 

momentum thickness 


viscosity 

^0 

permeability of free space 

P 

fluid density 

T e 

"edge" stress 

T w 

wall stress 

Subscripts: 


L 

laminar 

o 

entrance value 

T 

turbulent 
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w wall 

<=° asymptotic value 


ANALYSIS 

The problem concerns the flow of an incompressible, viscous, electrically conduct- 
ing fluid in the entrance region of a channel (see fig. 1). The fluid is assumed to have 
scalar conductivity and viscosity. The width of the duct in the Z -direction is taken to be 
infinite so that the problem is two-dimensional in the independent variables X and Y. 
Constant electric and magnetic fields are impressed in the Z and Y directions, re- 
spectively. The slip magnetic Reynolds number 

Re m = ^o a |f ' < U > 

is taken to be small so that the induced magnetic field is also small. Finally, the fluid 
is assumed to be electrically neutral and Hall effect-free. 

For the uniform entry case it is common to use the following boundary layer and free 
stream momentum integral equations (e. g. , refs. 10 to 12): 

Boundary layer: 


— (pU 2 0) + pU6* = t - ctB 2 U6* (1) 

dX dX w 


Free stream: 


= ctB(E - UB) - pU — (2) 

dX dX 

Maciulaitis and Loeffler use equations (1) and (2) even for the nonmagnetically fully de- 
veloped inlet condition. This is clearly a questionable procedure since, in this case, 
there is no reason to suspect that the pressure gradient is accurately given by equa- 
tion (2) near the inlet. Schlicting (ref. 7) and Roidt and Cess (ref. 9) also use equa- 
tion (2), but since theirs are Blasius series solutions, the unintegrated form of equa- 
tion (1) is used. 

Clearly, equations (1) and (2) cannot hold for large X because (1) the asymptotic 
solution t w = ctB 2 U6* is not, in general, compatible with the correct asymptotic solu- 
tion and (2) there is no shear stress term in equation (2). Both of these difficulties stem 
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from the same source. Equations (1) and (2) are only valid in regions where a true 
boundary layer exists (i. e. , a region where the Y-extent of the constant-velocity free 
stream is very large in comparison with the boundary layer thickness) . The region near 
the inlet of the uniform entry case is an example of such a region. Far from the inlet 
the boundary layer thickness and free stream width can be of the same order of magni- 
tude. The result is that a velocity profile and shear stresses develop across the entire 
channel in the real flow. In this way the downstream pressure gradient becomes depend- 
ent on the wall stress. If the boundary layer - free stream model is to produce a similar 
dependence, allowance must be made for a viscous coupling between the boundary layer 
and the free stream. This is done in the present work by permitting an edge stress to 
exist at the boundary layer edge Y = 6. 

The manner in which the edge stress enters can be seen by integrating the MHD ver- 
sion of Prandtl's boundary layer equations: 

pu ^L + pv^-= iL-dp + aB[E -uB] 

9X 9Y 3Y dX 

*L+ = o 

0X 0Y 

The integration is performed over the variable Y for both the boundary layer and 
the free stream. To compute the integrals of the inertia and Lorentz force terms in the 
free stream, the approximations u = U(X) and v = 0 are used. Also, a shear stress 
t e is allowed to act at the boundary layer edge. The integrated equations of motion are 
the following (e = 8/a): 

Boundary layer: 


— (pU 2 e) + pU6* = T — — - CTB 2 U6* (3) 

dX dX w (1 - e) 


Free stream: 

^ = <JB(E - UB) - pU — - - — - — (4) 

dX dX a(l - e) 

As noted, the derivation of equation (4) assumes the velocity U is constant over a cross 
section of the free stream flow. This may be bothersome to the reader in that the exist- 
ence of shear stresses (in particular, an edge stress) in a flow demands that the velocity 
profile be nonconstant. The velocity U in equation (4) is, therefore, more aptly thought 
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of as an average value of velocity over the cross section. However, for computational 
purposes, it is taken to be invariant with Y. 

The asymptotic condition is now 

fr - — — - tfB 2 US*| = 0 ( 

L (1 - e) Jx=°° 

Substituting equation (5) in equation (4) and noting that dU/dX approaches zero as X 
becomes large gives 



Equation (6) is the correct form for the balance of pressure, shear stress, and electro- 
magnetic forces in the fully developed flow. If, further, boundary layer - free stream 

model is chosen so that (1) the fully developed wall stress r is the correct value 

w, 

and (2) the mean velocity 



is the correct value, then the fully developed pressure gradient is automatically correct. 

It is clear that, in order that equations (3) and (4) be valid for all X, the edge 
stress r must vary. Since, at present, the only knowledge of t concerns its end 

C v 

values, assume momentarily that T g can be put in the form 



Equations (3) and (4) then become the following: 
Boundary layer: 


— (pU 2 0) + pU6* ^ = t 
dX dX 


w 


1 - f 


crB 2 U6* ’ 


w 


( 8 ) 
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Free stream: 


^2- = ctB(E - UB) - pU — - 
dX dX a 


( 9 ) 


The use of equations (8) and (9) now awaits only the selection of the boundary layer pro- 
files and the function f. This will be done for laminar and turbulent flow separately. 


Laminar Flow 

The laminar flow profile is assumed to be parabolic: 

* = 2 v -v 2 ( 10 ) 

U 


where tj - Y/S. 

This profile is chosen because it is simple and, for zero Hartmann number, exactly 
describes the true, downstream Poiseuille flow. 

The wall stress is then given by 


T =2pU 
w 6 


( 11 ) 


The derivative du/dY is zero at the boundary layer edge Y = 6. Although the idea of an 
edge stress r 0 at a point where the velocity derivative is zero may seem contradictory, 
it should be remembered that the boundary layer - free stream model is an approxima- 
tion to the true flow, which contains shear stresses throughout. Because of this approxi- 
mation the relation t = p (du/dY) cannot be used to calculate the viscous stresses at all 
points in the flow. In particular, it cannot be used at the boundary layer edge or in the 
free stream. The requirement of a nonzero stress at the boundary layer edge is not, 
therefore, a contradiction. It is a necessary condition if the boundary layer - free stream 
model is to yield accurate results both near and far from the entrance. 

When equation (10) is used, the quantities 9, 6*, and U can be computed in terms 
of e, a, and(U): 


9 = 




15 
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- <u) 

i-£ 


Inserting these equations together with equation (11) into equation (8) yields 


■ 2 + i 6 




d(e*) 2 _ 120 
dx Re 


2(1 - f) _ 


(e*)‘ 


( 12 ) 


where 



Re = 4 P< u > a 
M 



a 


At this point it is necessary to decide the manner is which the asymptotic boundary 
layer - free stream solution shall approximate the true Hartmann flow. Above all, the 
asymptotic wall stress and Hartmann wall stress should be equal. The latter is given 
by (ref. 13) 



Hartmann 


pMjj(U)tanh Mjj 



tanh MjA 
M H / 
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The asymptotic boundary layer - free stream wall stress is given by equation (11). 
Noting that 6 = e a and U = (U)/(l - e/3) results in 


2m(U) 


w, 00 




ae. 


Equating the previous quantities results in the quadratic equation 


tanh M n 


6 1 - 


H' 


M, 


C-3€ co + 


H 


M H tanh 


= 0 


(13) 


The appropriate root of equation (13) for the present problem is 



(14) 


The selection of the negative root is based on the fact that e must equal one when 
Mtt = 0. The value of f at X = °° is obtained by equating the right-hand side of equa- 

XI 

tion (12) to zero: 


M- 


H 2 


‘C 1 . ^ O * f « - 1 - -r 6 


(15) 


The initial values of f depend on the inlet conditions. The two obvious cases are 
the (1) uniform velocity and (2) fully developed (Poiseuille flow) profiles. The values of 
e and f at the inlet are, therefore, as follows: 

Uniform velocity inlet: 



(16) 
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Fully developed inlet: 


e 


o 


= 1 




( 17 ) 


In the uniform velocity case, f = 0 so that the flow near the inlet is Blasius-like. 
In the fully developed case, f is one for zero magnetic field strength (eq. (15) for 
m h = 0) . 

Equations (15) to (17) summarize all that is known of the function f. This author 
has not been able to deduce the functional relation of f to the other variables in the 
problem although this may, indeed, be possible. Instead, an inductive approach which 
uses certain known limiting solutions is used. The function f is required to be such 
that the solutions of equation (12) are consistent with these limiting solutions. The re- 
quirements are now listed: 

(a) The function f should be such that e varies monotonically with x. That is, if 
> e Q , e should monotonically increase, and, if e M < e Q , e should monotonically 

decrease. 

(b) For large x-values, e* should be expandable in the form (refs. 7 and 9) 


e* = 1 - a^exp(-Ajx) + . . . 


(18) 


(c) For the uniform inlet case at low x-values, e* should be expandable as (refs. 7 
and 9) 


- M 1/2 (i 


b 1 + 


V 1/2 ♦ 


°3 X 


+ . 


(19) 


when is of order 1. 

(d) For the uniform inlet case at low x-values, e* should be expandable as (ref. 4) 

e* = (x)*/^! + CgX + CgX^ + . . (20) 

when e w is very small compared to unity. 

(e) As x — °°, the value of e* should approach unity. 

For the case f = 0, the solutions of equation (12) satisfy conditions (a) to (d). It is, 
therefore, logical to choose f so that the right-hand side of equation (12) has the same 
kind of dependence on e* as it does when f = 0. However, when f = 0, the resulting 
solutions do not satisfy condition (e). One choice of f which has the appropriate depend- 
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ence on e*, but which also satisfies conditions (e), is obtained by demanding that the 
right-hand side of equation (12) be directly proportional to (1 - e*2) or 


120 

Re 


2(1 - f) _ % e *2 


= O'. 


(1 - e* 2 ) 


( 21 ) 


With this choice, the solutions of equation (12) asymptotically approach the condition 
e* = 1 when X becomes large. The proportionality constant is chosen by insisting 
that equation (21) be compatible with the entrance conditions (e. g. , eq. (16) or (17)). 

The result is 


a 


L 



2d - f 0 > 

-^i e *2 


2 

3 ° 

120 

^OO 


Re 

1 - 

e *2 



o 


If one uses this definition of a equation (21) can be put in the form 


f = f + 
o 


Af 

Ae*^ 


( e * 2 - e o 2 ) 


(22) 


where 


Ae* 2 = 1 - e* 2 

It is clear that a more arbitrary choice of f would not satisfy all of conditions 

(a) to (e). For instance, a choice which makes the right-hand side of equation (21) equal 

to at (1 - e*) would satisfy every condition but (d). This does not preclude other possi- 
-Li 

bilities, of course. However, the one given by equation (21) is reasonable, simple, and 
provides accurate analytical solutions. 

The criteria (a) to (e) say nothing of the near entrance behavior of the fully developed 
entrance case. This is because, to this author's knowledge, no such information exists. 
Therefore, the behavior of f according to equation (22) does not necessarily apply to 
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this case. However, since the chosen f is not necessarily inapplicable and since no 
other solutions are available, calculations are performed for this case also. 
Substituting equation (21) gives equation (12) as 



which can be rearranged to give 



Equation (24) can be integrated by partial fractions to yield 
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<P 3 = -cp 1 - cp 2 


* 4 = 



Using equation (11) and the continuity relation 
friction factor in the form 


U= (U)/(l - 6^6*73) gives the local 


2t. 


C fL " 


w 


16 


P< U) 


Re 1 - 



(26) 


The average friction factor ( C f is given by 



where 


*1 



^2 
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- *1 



Finally, the dimensionless pressure defect is obtained by solving for r /(I - e) in 
equation (3), substituting into equation (4) using the relations 6 = ea and U = (U) /(I - e/3), 
and integrating the resulting equation. This gives 



where K=E/(U>B. 

Turbulent Flow 

The entrance length problem for turbulent flow is complicated by the fact that no 

simple formula for the fully developed flow is known for arbitrary magnetic field 

strengths. Murgatroyd (ref. 14) provides experimental results in the Reynolds number 
4 5 

range 10 < Re < 10 . His experiments show that the wall shear stress varies from the 
Blasius formula value at zero Hartmann number to the Hartmann value at high Hartmann 
numbers. The behavior is such that, except for a small region near Mjj = 0, the wall 
stress monotonically increases with Hartmann number. Harris (ref. 15) has done an ex- 
cellent dimensional analysis of the fully developed flow and shows that the fully developed 
wall stress depends specifically on the parameters Re, Mjj/Re, and M^yfee. Unfortu- 
nately, he is unable to evaluate certain unknown functions of these parameters for all 
magnetic field strengths. Therefore, rather than attempt an analytical description of the 
downstream flow, the results of Murgatroyd (ref. 14) are used to obtain t w oo . 
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For present purposes it is assumed that the velocity profile and wall stress are 
given by 


— = tj 1/7 (29) 

U 


and 


T - fr>U 2 / M \ 1/4 

w 2 1/2 \pU5/ 


(30) 


where /3 = (0. 0225) 2 1//2 = 0. 0318. 

By using equations (29) and (30) in equation (8), one obtains, in a manner similar to 
that already outlined for laminar flow, the governing equation for turbulent flow: 



7, £ \ 5/4 1 

1 

d ( e *) 5 / 4 _ 90 

H) 

/3(1 - f) 4M H L _ e\ 3/4 r *5/4 

dx Re 

7 + 2e 

£ 5/4 Re 3/4 V 8/ 


Letting 

U ') 5 /4 . A _ i ^ 3/4 £ * 5 / 4 


using the approximations 


d(e’) 


5/4 


l 1 




3/4 


de* 


5/4 



7 + 2e 



(31) 


(32) 


fr 


16 



I 


* 


and assuming that 


gives equation (3 1) as 


where 


f-f + ** .. frS/4 _ e ,5/4\ 

° Ac’ 5 / 4 ' ° J 



(33) 


(34) 



The solution to equation (34) is 

e' 5 / 4 = 1 + ^q 5//4 - l^exp(-a' T x) (36) 

Unlike the laminar flow case the variation of the function f (eq. (33)) is not easily justi- 
fied. The only justifications offered here are that the variation given strongly resembles 
equation (22) for the laminar flow and the simple form of the solution given by equa- 
tion (36) would be impossible otherwise. Also, any variation of f that satisfies both end 
conditions is better than simply omitting f entirely, which is what is inadvertently done 
in references (5) and (11). 

As in laminar flow the two entrance conditions given by equations (16) and (17) apply 
for turbulent flow as well. The value of is obtained by equating the known value of 
T to the value given by equation (30) . The local friction factor C fT , is approximated 
by 
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with e approximated from equation (32) as 


e = 



(38) 


Equation (38) is slightly inaccurate at small x-values but approaches the value given 
by equation (32) as x increases. 

The average value (C fT ) is obtained from 
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where 
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Finally, the pressure defect is given by 



DISCUSSION 
Laminar Flow 

Figures 2(a) and (b) show the laminar flow velocity developments for = 0 and 
Mjj = 4 with uniform entrance conditions. These are compared with the numerical solu- 
tions of Bodoia and Osterle (ref. 6) and Shohet (ref. 8) . Even though the free stream 
velocities (large y-values) are somewhat in error, the velocities near the wall (small 
y-values) follow the numerical solutions quite closely. For friction factor calculations 



Figure 2. - Laminar flow velocity development. 
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only the accuracy of the velocities near the wall is relevant. Figure 2(b) also shows the 
development for the fully developed entrance case (Poiseuille flow) for = 4. The 
flow is as one might expect - deceleration at high y-values, acceleration at low y-values. 
No comparisons with numerical solutions are possible here since, to this author’s know- 
ledge, no such solutions have been obtained. 

Figures 3(a) and (b) illustrate the pressure defect development corresponding to 
figures 2(a) and (b). In figure 3(b), the load parameter was also varied and its effect is 
evident. The agreement is excellent. 

Figure 4 shows the hydromagnetic pressure defect development for = 4 
and both entrance conditions. The hydromagnetic pressure defect is defined as 

4M H 

+ x(l - K) (42) 

Re 

From equation (28) it can be seen that is independent of the electric field. Figure 4 
indicates a less severe entrance pressure drop for the fully developed entrance condi- 
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tion. It also shows that, for this case, a greater distance from the entrance is required 
for the pressure gradient to approach its asymptotic value. 

Figure 5 is a comparison of the present work with the numerical flat plate results of 
Dix (ref. 3) and the momentum integral - Hartmann profile method of Heywood and 
Moffatt(ref. 16). The solution given by equation (26) is representative of flat plate results 
only when the Hartmann number is very large. Otherwise, the free stream veloc- 
ity varies with X. The value M H = 100 was chosen as being representative of a large 
Hartmann number. It can be seen that the agreement is again quite good. 



Figure 6. - Entrance lengths for laminar and turbulent flows. 
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Entrance lengths for both entrance conditions are depicted in figure 6. The en- 
trance length X 0 is defined here as the distance from the entrance at which the corre- 
sponding e(x) satisfies 


e - e 0 = 0. 99(6„ - g (43) 

It is evident that the entrance lengths for the fully developed entrance case are siz- 
ably larger than for the uniform entrance case. The behavior is similar, however; both 
approach constant values as Mjj becomes small. 


Turbulent Flow 


When the Hartmann number is zero, equations (29) and (30) should be adequate to 
describe the turbulent flow entrance problem. This problem was examined (for the uni- 
form inlet case only, of course) and the following quantities were obtained: 

X e = 12. 1 Re 1 / 4 a (44) 

A^ t = 0. 088 (45) 


where X g is the entrance length, and is defined as 


A^p = lim 

x—°° 


/8 \7/4 / i \l/4 I 

(-) 4 —) - 

\7/ \Ree/ 1 


(46) 


The quantity A#p is therefore the difference between the dimensionless pressure drop 
in a very long section of channel where the flow is fully developed (first term in brackets) 
and the drop in an equal-length channel with entrance effects #p. 

Equations (29) and (30) were found to be inadequate for describing the flow at 
Hartmann numbers other than zero. Asymptotic wall shear stresses computed from 
Murgatroyd's (ref. 14) data frequently yielded negative edge stresses and, for some 
Hartmann numbers, boundary layer thicknesses greater than the channel half width. 

Both of these results are physically meaningless. 

The only reason for trying equations (29) and (30) at Mjj ^ 0 is simply expedience; 
no better models are available. However, the erroneous results mentioned in the pre- 
ceding paragraph indicate than any computations based on these results would also be 
erroneous. Hence, no calculations were performed. Clearly, if a meaningful descrip- 
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tion of the asymptotic flow cannot be found, there is no point in making entrance flow 
calculations. It should be apparent, however, that this is not a shortcoming of the edge 
stress method. 


CONCLUDING REMARKS 

A momentum integral method has been developed for describing the incompressible 
flow of a conducting fluid in the entrance of an MHD channel. The method is amenable to 
various initial conditions and the wall shear stress asymptotically approaches the correct 
downstream value. All this is possible by introducing a boundary layer edge stress, the 
correct choice of which forces the governing differential equation to the correct asymp- 
totic solution. The technique is accurate for describing laminar entrance flows, with or 
without a magnetic field. Turbulent flow without a magnetic field can also be analyzed. 

The analysis of turbulent entrance flow in a magnetic field is dependent on a know- 
ledgable choice of boundary layer velocity profiles and wall stresses. In the present 
work one-seventh power velocity profiles and Blasius formula wall stresses are used. 

This description is found to be good for a Hartmann number of zero but inadequate other- 
wise. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 7, 1969, 

129-02-01-13-22. 
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